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1. Introduction 

Despite contributing only a small fraction of the overall cosmic-ray flux, cosmic-ray electrons 
and positrons (CREs) provide an important and unique probe of our local Galactic neighborhood. 
CREs lose energy rapidly via inverse Compton scattering and synchrotron processes while prop¬ 
agating in the Galaxy, effectively placing a maximal propagation distance for TeV electrons of 
order ~1 kpc[[Ij]. CREs at TeV energies provide a direct measurement of the local cosmic-ray 
accelerators and diffusion. 

The Fermi-LAr[J|], and more recently A MS [[3]]. have measured the CRE spectrum with high 
statistics using satellite-based experiments up to energies of several hundred GeV. However, above 
those energies these instruments run into difficulty from the combination of the steep CRE spec¬ 
trum and their relatively small acceptance. Ground-based experiments that utilize the Imaging 
Atmospheric Cherenkov Technique (IACT) have the capability to extend the CRE spectrum to 
higher energies due to their large collection area ( 10 5 m 2 ). HESS[||[|] and MAGIC [0] have 

demonstrated this ability and previously measured the CRE spectrum from the ground up to TeV 
energies. Their results provide evidence of at least one nearby CRE source and agree with satel¬ 
lite measurements within systematical uncertainties where there is energy overlap. The combined 
picture that has emerged is one where the CRE spectrum is mostly flat and can be described by 
a simple powerlaw from ~10 GeV up to just below ~1 TeV, above which HESS has measured a 
spectral steepening. MAGIC data[]6j] are consistent with a single power-law up to ~3 TeV. 

Additionally, the positron fraction spectrum, 0(e + )/(0(e _ ) + 0(e + )), has now been measured 
above 10 GeV by the HEAT|j7|], PAMELA[|J], Fcnni-L47’[fi]], and AMS[|Io|] experiments. This ra¬ 
tio is found to rise with increasing energy up to ~200 GeV, above which it appears to flatten out. 
Positrons are mainly produced in secondary interactions between cosmic rays and the interstellar 
gas. These results point to the possible existence of an additional local source of positrons on top of 
secondary production. Explanations could include additional production from nearby standard as- 
trophysical objects, such as pulsars or supernova remnants, or more exotic production mechanisms, 
such as the annihilation or decay of particle Dark Matter. Conversely, recent studies have also pro¬ 
posed a more detailed propagation model Jl2|] or a better accounting of secondary production[|T3]] 
to describe the results. A full understanding of this situation will require detailed input about both 
the positron fraction and the CRE spectrum. While the excitement of the unexpected excess found 
in earlier ATIC data[|TT| | is now largely over, a high statistics measurement of the CRE spectrum to 
TeV energies will help us build a clear picture of our local CRE emitters. 

2. Methods 

Data presented here were collected by the VERITAS telescope array located at the Fred 
Lawrence Whipple Observatory (FLWO) in southern Arizona (31° 40'N, 110° 57'W, 1.3km a.s.l.). 
VERITAS is a ground-based array of four telescopes sensitive to gamma- and cosmic-rays above 
~100 GeV. While VERITAS is primarily a gamma-ray instrument, CREs are a diffuse source 
across the sky and so are collected during all astrophysical observations. 

To isolate this signal, we apply strict analysis cuts to select only the best reconstructed events 
under pristine weather conditions. Data events must have four good telescope images and recon- 
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BDT Response 


Figure 1: BDT response for proton MC (gray filled area) and the full dataset (blue points) for the energy 
range 630 GeV to 1 TeV. The insert shows the ratio, data/MC, over the same range. The agreement is very 
good except close to the limits of the distribution. As we approach the side dominated by background-like 
events, -1.0, we find an excess in the data over the proton MC. We expect an excess here from helium and 
higher-Z primaries, particularly since helium makes up ~ 20% of the cosmic-ray flux. Likewise, as we 
approach the side dominated by signal-like events, 1.0, we find an excess in the data over the proton MC that 
arises from the CREs measured in this study. 


struct within the inner 1° of the camera (the VERITAS field of view is 3.5°). A strict cut is placed 
on the distance between the array center and the reconstructed array core, coreR < 200 m. Only 
extragalactic fields are considered in order to reduce the contamination from the Galactic Plane 
diffuse gamma-ray flux. Additionally, all detected or candidate gamma-ray sources within each 
field of view are excluded. To reduce the number of detector configurations in our simulations, 
we selected only data collected between September, 2009 and July, 2012, which are the dates of 
the two major VERITAS hardware upgrades (telescope relocation and PMT replacement, respec¬ 
tively). We also restrict the mean zenith angle of the data to be between 65° — 75°, where we 
exclude data runs with mean values outside of this range (a data run is typically ~20 minutes). To 
accumulate sufficient Monte Carlo statistics, we generated simulations at a single zenith angle, 70°. 
This restricted data zenith range ensures the level of data/MC agreement necessary in this analysis. 
296 hours of live-time remain after all these cuts. 

We rely heavily on our Monte Carlo for interpretation and signal extraction. To simulate 
the electromagnetic and hadronic showers, we used Corsika 6.970[|l6|], with the QGSJetII.3 and 
URQMD 1.3cr underlying event generators, and GrISUDet 5.0.0[|T7|l for the VERITAS detector 
response. We generated electrons, protons, and helium showers with a 4° radius on the sky to 
approximate the isotropic and diffuse cosmic-ray flux. Larger simulation radii were tested and 
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Figure 2: The top three panels present the fit results for each of the 15 trials in three energy bins. The x-axis 
represents the trial number and the y-axis represents the measured electron fraction (the individual trial error 
bars are the let errors from the fit). Also shown is the mean value and uncertainty on the mean for each bin. 
The bottom five panels show the individual fit results for the first five trials in the 631 —708 GeV energy bin 
(the black points are data, the best fit line is wide and red, the proton MC line is thin and green, and the 
electron MC line is dashed and orange). For each fit we list the best fit electron and proton fractions and the 
X 2 /dof for the comparison between the data and best fit. 


found to not improve this approximation. 

For signal and background discrimination we use Boosted Decision Trees (BDTs) that were 
integrated into the standard VERITAS analysis chain using the ROOT TMVA[|T4|| framework. The 
BDTs were trained with a diffuse electron MC sample (signal) and a representative sub-sample 
of the data chosen randomly from the full dataset (background). We used four array-level shower 
variables in the training and event discrimination: MSCW, MSCL, ^ 2 (£), and the emission height. 
MSCW and MSCL are variables based on the comparison of the spread along the major and minor 
axes of an ellipse fit to the camera images (length and width respectively) compared to expected 
values from simulations. % 2 (E) represents the variability of the energy measurements in each of the 
four telescopes. The emission height is the reconstructed height of the peak width of the shower. 
We ensured that the trees were not overtrained by using a fraction of the full training signal and 
background samples to test the response. No overtraining was found and the trees and method were 
also used to successfully reconstruct the Crab Nebula energy spectrum as a cross-check. 

Each data event is assigned a BDT response value, from — 1.0 to 1.0, where higher values in¬ 
dicate that the event is more signal-like. Figure [I] shows a comparison of the BDT response for the 
full dataset and proton MC. The agreement is very good except near the limits of the distribution. 
As we approach the side dominated by background-like events, —1.0, we find an excess in the data 
over the proton MC. We expect an excess here from helium and higher-Z primaries, particularly 
since helium makes up ~ 20% of the overall cosmic-ray flux[]T5[]. We investigated the BDT re- 
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sponse of helium MC events and they are found to peak at —1.0 and fall off faster than proton MC, 
in agreement with this interpretation. As we approach the side dominated by signal-like events, 
1.0, we find an excess in the data over the proton MC that arises from the CREs measured in this 
study. 

We select only those events with BDT response values > 0.7, which focuses on the region 
that contains the majority of the signal-like events and rejects the majority of the background-like 
events. We then employ an extended likelihood fitting method to the BDT response within this 
region to extract the contribution of electron and proton events to the total. This fit floats the 
electron and proton MC shapes relative to each other to find the best combined fit to the data. 
Helium and higher-Z shower events are found to be sufficiently rejected by the BDT cut to ignore 
at first order. To estimate the final electron fraction and its uncertainty, we divide the data into 
sub-samples and run the experiment several times (see Fig. |2|). The mean of these separate trials 
is the final electron fraction with the uncertainty defined as the error on this mean. This method 
removes some systematic biases from the fitting method and lets the data itself drive the statistical 
precision of the measured value. One caveat of this technique is that fits in trials/bins with less than 
~100 data events do not return sensible results. As a result we use fewer trials for higher energy 
bins where there are less statistics. The highest bin presented here is the result of a single fit to 200 
events; the uncertainty quoted is the la error from the fit. 

3. Results 

We show in Fig. |] the preliminary VERITAS CRE energy spectrum spanning ~300 GeV to 
~5 TeV. The spectrum steepens at higher energies and is best described by two power-laws with a 
cutoff. The best fit for this cutoff energy is found to be 710±40 v/a/ GeV, with best fit spectral indices 
below (above) this energy of — 3.2 ± 0.1 i/a , (—4.1 ±0.1 ifaf ). The X 2 /d°f of this fit is 9.71/11. The 
gray band represents the systematical uncertainty, which is dominated by the ~ 20% uncertainty 
on the VERITAS absolute energy scale. This translates into a +64%/—33% (+98%/—43%) system¬ 
atical uncertainty for a spectral index of —3.2 (—4.1). We measure an additional 10% systematical 
uncertainty above ~ 1 TeV that quantifies hardware uncertainties at those energies. This additional 
uncertainty is added in quadrature. 

Due to the similarity of gamma and electron electromagnetic showers, we cannot rule out a 
significant contamination of gamma-ray events within our electron data. However, Fermi-LAT has 
now measured the diffuse extragalactic gamma-ray flux up to ~800 GeV and finds it orders of 
magnitude below their own CRE measurement!^ ]. They additionally find evidence for a cutoff in 
the diffuse gamma-ray spectrum above a couple hundred GeV. 

CRE results shown here qualitatively agree with prior ground-based and satellite-based mea¬ 
surements at similar energies. Of the many experiments studying CREs, this result represents the 
second high-statistics measurement of a cutoff in the CRE spectrum around ~ 1 TeV. The precise 
measurement of this energy cutoff is an important parameter in any successful model our local 
CRE environment. VERITAS has significantly more data on disk than what has been used in this 
study and work continues with the goal of extending our measurement out to even higher energies. 
The CRE spectrum between 5 and 10 TeV is unexplored and this is the energy range where we ex¬ 
pect to see spectral features from individual nearby astrophysical sources (if they are the dominant 
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Figure 3: VERITAS preliminary cosmic-ray electron spectrum as a function of energy (GeV) in solid blue 
circles. The best fit to the data is represented as a dashed line and is found to be two power-laws with a break 
energy 710 ± 40 GeV and spectral indices of —3.1 ± 0.1 stat (—4.1 ± 0.1 stat ) below (above) the break. The 
^ 2 /dof of the fit is 0.9. Shown for comparison are data from other experiments in the same energy range. 
The gray band represents the systematical uncertainty on the VERITAS measurement. 


source). We urge caution in over-interpretation of the uptick in the final VERITAS data point since 
this is within 2a of the best fit line. 
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